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ABSTRACT 

We investigate the dust size and dust shell structure of the bipolar proto-planetary nebula M 1-92 by means of radiative transfer 
modeling. Our models consists of a disk and bipolar lobes that are surrounded by an AGB shell, each component having different dust 
characteristics. The upper limit of the grain size a max in the lobes is estimated to be 0.5 yt/m from the polarization value in the bipolar 
lobe. The a max value of the disk is constrained with the disk mass (0.2 M©), which was estimated from a previous CO emission line 
observation. We find a good model with a max = 1000.0 /um, which provides an approximated disk mass of 0.15 M©. Even taking into 
account uncertainties such as the gas-to-dust mass ratio, a significantly larger dust of a max > 100.0 //m, comparing to the dust in the 
lobe, is expected. We also estimated the disk inner radius, the disk outer radius, and the envelope mass to be 30 i?*(=9 AU), 4500 AU, 
and 4 M©, respectively, where v exp is the expansion velocity. If the dust existing in the lobes in large separations from the central star 
undergoes little dust processing, the dust sizes preserves the ones in the dust formation. Submicron- sized grains are found in many 
objects besides M 1-92, suggesting that the size does not depend much on the object properties, such as initial mass of the central star 
and chemical composition of the stellar system. On the other hand, the grain sizes in the disk do. Evidence of large grains has been 
reported in many bipolar PPNs, including M 1-92. This result suggests that disks play an important role in grain growth. 

Key words. Stars: AGB and post- AGB - circumstellar matter - radiative transfer - polarization - individual (M 1-92) 



1. Introduction 

Planetary Nebula (PN) morphology represents a history of 
physical processes in the stellar/ circumstellar environments. 
A large sample (-500) of galactic PNs has shown that 
20-25 % of them have spherical structures, but the oth- 
ers have asymm etric (elliptic, bipo l ar, an d point- symmetric 
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Several theoretical studies have 
been car ried out to exp l ain th e PN morphologies (see re- 
view by IBalick & Fra nk 2002). Widely accepted facts and 
physical proc e sses are cool spots o n the stellar surface 
iFrankl 119951: ISoker & Clavtonl 119991) . the ste llar rotation 

magnetic field 
1998, [19991: 



(Biorkm an & Cassinellil 11 993). 
19921: iMastrodemos & Morris! 



Pascoli et al. 
Matt et al.l 12000). 



and binary interaction (iLivio. Salzman. & Shaviyll 19791: iMorrisI 
I1981L 119871: iBond & Livid 119901: ISoker & Livid 11994 . While 
details are still unknown, these mechanisms one way or an- 
other cause the equatorially enhanced mass loss in the AGB 
phase (e.g. iMeixner et all 1 19991: lUeta et al.l 12000). In the sub- 
sequent post- AGB phase, the central star blows away the low 
density fast wind, which inflates the l obe in the polar direction 
more than in the equatorial direction (iKwok et al .11 197 8t iKwokl 
Il982l) . This interacting stellar wind (ISW) mo del explains y ari- 
ous morphologies of elliptic and bipolar PNs (lBalicklll987l) . In 
the above scenario, the binary interaction is thought to be one 
of the most promising mechanism to shape PNs with a narrow 
waist between the bipolar lobes because this can amplify the 
mass loss in the equatorial plan e and make a disk-like structure 
(e.g. ISoker & Rappaport 2000). In addition, the binary interac- 



tion is also thought to form a symmetric structures such as jets, 
ansae (e.g.lSoker 1990. 119921) . and point- symm etric shapes such 
as a spiral and quadruple by pro cessing motion ( Manc hado et al.l 
ll996UMauron & Huginsll2006l) . 



It is difficult to directly detect the effect of binary interaction 
on PN shaping using the current observing techniques. However, 
it would be possible to infer the presence of the effects due to the 
aforementioned scenario if the inner part of the circumstellar en- 
vironment is investigated closely enough. Near-infrared imaging 
polarimetry is a powerful method to probe dust shells and to pro- 
vide important information that can be derived in more details 
than non-polarimetric imaging. Previous observations and dust 
scattering modeling have re vealed disk and bipolar lobe struc- 
tures of evolved stars (e.g. IScarrott"&S carrott 1995; ISu et al.l 
2003; Gledhi O 120051: lUeta et al.l 120051: iMurakawa et al] 120051 
l2007[lUeta et al.ll2007[ hereafter UMM07). However, most pre- 
vious experiments treated only a single grain model in the en- 
tire dust shell or SEDs were not considered simultaneously. 
These shortcomings have prohibited us from investigating into 
the grain and disk properties in detail. To resolve the situ- 
ation, we considered both polarimetric images and SEDs in 
our recent radiative transfer modeling of the bipolar proto-PN 
(PPN) Frosty Leo, assuming two different dust models, one in 
the equatorial region (d isk) and the other in the bipolar lobes 
(Mura kawa et al.ll2008al) . From this modeling, we were able to 
derive evidence for grain growth in the disk region. In this pa- 
per, we use this technique for the bipolar PPN: M 1-92 (=IRAS 
19343+2926), whose characteristics were discussed previously 
by us (UMM07). In Sect. [2 details of our radiative transfer mod- 
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eling are described. We will discuss the implication of our model 
result in Sect. 

2. Radiative transfer calculations 

The primary purpose of our modeling is to explore the parameter 
space of the disk geometry and the grain sizes in the circumstel- 
lar dust shells of our target. In order to do this, we applied model 
geometries with a disk and bipolar lobes and different dust mod- 
els for them. We used our own three- dimensional Monte Carlo 
code STSH (Mu rakawa et al.l l2008bl) . which solves radiative 
transfer problem of scattering and absorption by dust. This code 
can handle multiple dust models in an arbitrary model geometry 
and computes SEDs, dust temperatures, and the Stokes IQUV 
parameters. We constrained the model parameters by compar- 
ing with the observed SEDs and polarimetric images from the 
HST NICMOS 2 data at a wavelength of 2 u rn- We followed 
a simi lar modeling procedure as in the past (iMurakawa et al.l 
2008a : b). We first tried several parameter sets to find an approx- 
imated solution and parameter ranges and determined some pa- 
rameters, which can be estimated easily. The grain sizes in the 
bipolar lobe were determined from the polarization values in this 
region, which do not depend on the other parameters such as 
the model geometry and the masses of the disk and envelope. 
Then, we constructed a large number of models by SED fit and 
made model images of a few possible parameter sets. We find 
some good models with different dust sizes a max in the disk, 
which reproduces the characteristics of the above observations 
well. To constrain the <2 max value, we use the disk mass which 
was estimated from a previous CO emissio n line observation 
dBuiarrabal et al]ll998tt lAlcolea et aTll2007l) . We will describe 
this grain size issue in detail below. 



Table 1. Model parameters of our radiative transfer calculations. 



parameters adopted values comments 





central star 






6500 K 


2 


U 


7000 L Q 


adopted 3 

4 


d 


3.0 kpc 


R* 


1.3 x 10 12 cm 


calculated 




disk 




Rm 


30Z?*(=9AU) 


30-50 


Rdhk 


4500 AU 


3000 - 4500 


H 


0.3 


0.3 - 0.4 


T2.0 


30 


adopted 3 


M disk 


0.15 Mq 


5 


^max 


1000.0 fim 


> 10.0 


./in 


0.2 


adopted 3 




superwind and AGB shell 




7/5 


adopted 6 


a 


1.2 


adopted 6 
adopted 6 




1.7 


7 


0.7 


adopted 6 


e m 


0.01 


adopted 6 


^rim 


0.5 


adopted 6 


^out 


30" 


30-60 


^max 


0.5 jjm 


adopted 7 


fm 


0.2 


adopted 3 


-^env 


4M Q 


4-6 


M 


7.5xl0- 6 [v exp kms" 1 ] M Q yr~ hS 



1 Rang es give the uncerta inty of the corresponding model parame- 
ters. 2 [Arrieta et al] ( |2005|). 3 Based on com parison of the SED, 4 
ICohen & Kuhil dl977MEiroa & Hodappl d!989h . 5 with assumption of a 
^max = 1000.0 iim grain model, 6 Based on comparison of the intensity 
image, 7 Based on comparison of the polarization image, 8 calculated. 



2.1. Model assumption and numerical simulation 

M 1-92 is an oxyg en-rich PPN known as Minkowski's footprint 
( Mink owskil 1 9461) An optical spectrum suggests a 6500 K black 
body component and an 1 8 000 K component in shorter wave- 
leng ths. The latter is p robably due to flux from the companion 
star (lArrieta et al .1120051) . Because the 6500 K component dom- 
inates in the total flux, we assumed a blackbody spectrum with 
this temperature as the illumination s ource. The distance D was 
determin ed to be between 2.5 kpc (ICohen & Kuhil [T977I) and 
3.5 kpc (lEiroa&Hodapplll989l) . We adopted 3 kpc, the av- 
erage of the two. For the luminosity, although ICohen & K uhi 
(1 19771) obtained a 1600 [D kpc] 2 L©, which corresponds to 
-20000 L Q at 3.5 kpc, we found that 7000 L fits the absolute 
flux in the SED better and adopted this value in our modeling. 

The nebula has a w aterlily petal-like shape with an 11 " x 
6" e xtension (see also Tram mel! et al.l 1 19961: iBujarrabal et al.l 
1998a), which extends towards the northwest and the southeast 
directions. The northwestern lobe is brighter than the southeast- 
ern one. Blue- and red- shifted components are detected in the 
northwestern and southeastern lobes, respectively, in the CO 
emission l ine data and the northw estern lobe is tilted towards 
us by 35° (Bujarrab al et al.| [l998b). In our modeling, we assume 
that the model geometry consisting of three components: (1) a 
disk in the innermost part, (2) a superwind shell in the outside 
of the disk, and (3) an AGB dust shell i n the outermost part. 
The disk form chos en is one studied by iToomord ([1982) and 
iFischer et all (119961) . For the superwind shell, the NICMOS 2 
intensity image shows that the surface brightness at the outline 
of the bipolar lob es is slightly enh anced compared to the pro- 
jected inner part (lUeta et al.ll2007l) . suggesting a hollow struc- 



ture, which is probably formed by the interaction with the ISW. 
The AGB shell is assumed to have a one-dimensional, spheri- 
cally symmetric structure with a single power-law radial density 
gradient. The mass density distribution p consists of the individ- 
ual components of the disk pdisk, the superwind shell p sw , and 
the AGB shell Pagb and is given by 



P - Pdisk + Psw + Pagb, 



PdiskO,z) = PdO/#in) exp 



Hr 



for R\ n <R< R< 



disk? 



Psw 6) - Pe (*/* 8w ) | for yRiQhQ < R < R ^ 

Pagb (R) = Pe (R/Rsw)~ 2 for R sw <R< R ouU 

where the coordinates of (R, 6) and (r, z) = (R sin 6, R cos 6) are 
the two-dimensional spherical coordinate and cylindrical coor- 
dinate, respectively. We assume that the disk exists in a region 
between the inner radius R m and the outer radius ^disk- The value 
H is the ratio of the disk height to the disk radius. The pa is the 
density coefficient and is derived from the optical depth in the 
equatorial plane instead of the disk mass M^ sk as done before 
dMeixner et al.l 120021: IMurakawa etaD 12008 ah . We find that an 
optical depth of 30 at a wavelength of 2.0 /im fits well in the 
optical to NIR fluxes. We apply this value for further investiga- 
tion of the grain size effect. The shape of the superwind is deter- 
mined by R^q = R sw QO/n - l/2\f ( Oppenheimer et al 1 l200l 
and references therein). The density factors of £i n and e^m pro- 
duce the hollow structure of the lobe. In the above parameters, 
we adopted an a of 1.2,/3 of 1.7, y of 0.7, ei n of 0.01, e r i m of 0.5, 
and R sw of 7 r /5, which fit well the appearance of the bipolar lobe. 
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position (stellar radius) position (stellar radius) 



water ice silicate water ice 




Wavelength (yLtm) Wavelength (/im) 

Fig. 1. Mass density distribution of the selected model (top left panel), dust temperature distribution (top right panel), dust opacities 
as function of wavelength (bottom left panel), and comparisons of model SEDs with observations (bottom right panel). For the 
selected model, the attenuated stellar flux, thermal emission, and scattered light from the dust are also plotte d with the dashed-curve, 
the dashed-dotted curve, and the dotted curve, respectively. The observed data are from Mermilliod (1991) and USNO catalog for 
optical photometry, point source catalogs of 2MASS, MSX, and IRAS for infrared photometry, and ISO SWS and LWS for mid- 
and far-infrared spectra, respectively. In the opacity plot, the mean opacities of absorption k, scattering cr, and extinction k + cr are 
indicated with dashed, dotted, and solid lines, respectively. Models with four different a max values of 0.5 /mi, 1.0 /mi, 10.0 yum, 
and 1000.0 jim are compared. The a max = 0.5 jim model, which is used in the superwind and AGB shells, is also presented for 
comparison. In the SED plots, several curves with different colors denote model results with different grain sizes (a max ) in the disk. 



The AGB shell exists in the outermost part between the radii of 
R sw and R 0VLt . The density coefficient p e is determined with the 
envelope mass M env . We applied a gas-to-dust ma ss ratio of 160 , 
which is often used for oxygen-rich evolved stars (Knapp ll 19851) . 
The model geometry is viewed at a viewing angle # v i ew measured 
from pole-on. The free parameters and their ranges are an R[ n of 
30, 50, and 100 R±; H of 0.2, 0.3, and 0.4; M env of 2, 4, 6, and 
8 M Q ; and <9 view of 65°, 55°, and 45°. 

For the dust grains, we simplify the model as much as pos- 
sible to focus on determining the grain sizes which character- 
ize the optical properties of the dust the most. In our model- 
ing, we assume spherical cores with a mantle and an MRN-like 



size distribution of 0.005 /im < a < a max jim with n(d) oc a 3,5 
(Mathis, Rumpl & Nordsieck 1977). The ^max values are differ- 
ent for the disk and the envelope. The chemical co mpositions of 
the core and the mantle ar e astronomical sili cate (iDrainell 19851) 
and crystalline water ice (iBertie et al.lll969l) . respectively. The 
water ice mantle is assumed to have a constant thickness pro- 
portional to the core radius, which is determined with a volume 
fraction of the water ice to the grain core f m . This value is chosen 
to be 0.2 which fits the 3 jim water ice absorption feature. The 
grain size in the lobe was estimated from the polarization values 
in the lobe, because the NIR polarization strongly depends on 
the grain size in cases of submicron size. In our NICMOS2 data, 
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Fig. 2. Model results of the intensity (upper panels) and polarization (bottom panels) images of Ml-92 comparing with the 
HST/NICMOS 2 data from UMM07 (left column). Models with grain sizes in the disk of d m2LX = 1*0 y^m (middle column) and 
^max = 1000.0 jim (right column) are presented to show the effects of grain sizes. In the polarization images, the polarization vector 
lines are also plotted. 



the 2 yum polarization is P = 40 - 50 % in the upper lobe and 
P = 50-55 % in the lower lobe. With a large inclination angle 
of 35°, which is expected from CO observations, the calculated 
polarization becomes too low in the upper lobe and too high in 
the lower lobe. The average polarization between the upper and 
lower lobes for the a max of 0.3, 0.5, and 0.7 yum are 60 %, 45 %, 
and 20 %, respectively. We apply the a max = 0.5 yum dust model 
in the envelope in the subsequent models. The free parameter 
is the a max of the disk and the values to examine are 0.5 yum, 
1.0 yum, 10.0 yum, 100.0 yum, 1000.0 /zm, and 10000.0 yum. 

2.2. Result 

From the aforementioned parameter sets, we selected a good 
model with the following parameters: R{ n = 30 H = 0.3, 
R disk = 4500 AU, M env = 4 M Q , R out = 30", # view = 55°, and 
^max = 1000.0 yum. These are summarized in Table [T] 

The top panels of Fig. [T] show the cross section of the mass 
density distribution and the dust temperature distribution. As 
seen in the mass density distribution map, the superwind shell 
has a waterlily petal- shaped appearance. In the dust temperature 
distribution map, intermediate to high temperature components 



(7d £ 100 K) are seen in the innermost low-density polar region. 
The rim of the superwind and AGB lobes has a low temperature 
of ^ 70 K, which characterizes the flux in the far-infrared (FIR) 
or longer wavelengths. 

The bottom panels of Fig.Q] show the opacity (i.e. the size- 
averaged cross section per particle mass) of the modeled dust 
and the SEDs. To show the effect of grain sizes, we also present 
other dust models with a max = 1.0 Jim and a max = 10.0 jjm for 
comparison. In the opacity plot, the difference of the wavelength 
dependence for different dust models is clearly visible. As a mSiX 
increases, (1) the wavelength where the opacity drops becomes 
longer and (2) the strengths of water ice and silicate features 
become weaker. The SED plot compares the model results of 
different grain sizes of a max = 1.0 yum, 10.0 yum, and 1000.0 jim 
to the observations. The model SEDs reproduce the water ice 
features at 3 yum, 44 jim, and 62 yum and the silicate features at 
9.7 yum and 18 yum. Because the disk temperature is higher than 
the temperature of the superwind and AGB shells, in principle, 
the grain-size effect appears in the mid-infrared (MIR) flux, and 
it is possible to constrain the grain size by SED fit. On the longer 
wavelength side, the difference in the SED appears only for A ^ 
50 jjm. 
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FigureO shows the model results of the intensity images (top 
panels) and polarization images (bottom panels). The results of 
two different dust models of a m2LX = 1.0 /mi (middle column) and 
^max = 1000.0 yum (right column) are compared with the obser- 
vation (left column). The intensity images reproduce a waterlily 
petal- shaped bipolar appearance reasonably well. In the polar- 
ization images, the difference of the polarization values between 
the upper and lower lobes is too strong, as mentioned before. 
With respect to the effect of the grain size, we do not find any 
sufficient difference from our results. 

Hence, it is difficult to determine the grain size in the disk 
of M 1-92 from the SED, intensity image, and polarization im- 
age, which we used in our modeling. The reason is explained as 
follows. The MIR flux F from the disk is approximately propor- 
tional to MdiskK (a)- Based on our assumption of a constant opti- 
cal depth of the disk for different grain models, the optical depth 
r is proportional to M^(k(o) + <x(<z)). Therefore, we obtain 
F octk (a) I (k (a) + o~ (a)). For grain models with a max > 10 fim, 
the k {a) I (k (a) + cr (a)) value becomes nearly constant. With the 
similar reason, the wavelength dependence of the dust scattering 
matrix elements become low in large grains, resulting that the 
effect of the grain size does not appear much in polarization im- 
ages in the optical to NIR. Thus, the presence of large grains 
are examined often using a spectral index of dust opacity ft in 
longer wavelengths, i.e. the su bmillimeter and millimeter wave- 
length ranges (e.g. I Jura et al .11 19971) . In objects with large grains 
in the disks, the spectral slopes of the flux in the se wavelengt h 
ranges are shallower due to lower values (e.g. lDraindl2006l) . 
However, since F and r depend on the disk mass, if the disk 
mass is determined by other mea ns, the grain size can be con- 
strained. iB^OT^aTetal] (11998b) obtained ~1" resolution im- 
ages of M 1-92 in the 13 CO 7=2-1 emiss ion line and an esti - 
mate of a gas disk mass of 0.2 Mq (see also lAlcolea et al .1120071) . 
The calculated disk masses in our modeling are 0.022, 0.0098, 
0.017, 0.049, 0.15, and 0.47 M Q for dust models of a max = 0.5, 
1.0, 10.0, 100.0, 1000.0, and 10000.0 /mi, respectively. We find 
that the a max = 1000.0 yum model provides the best estimate. 
However, facts such as the gas-to-dust mass ratio, a region that is 
considered to be the disk, and the dust size distribution can lead 
some uncertainties in the total disk mass. For example, the gas- 
to-dust mass ratio, which we adopted in our model is es timated 
from a number of mass-losing oxygen-rich giant stars ([Knapp 
119851) . The value o f individual objects v aries with a factor of ~2 
in general (see also lOlofsson et al.ll 19931) . In a special case of ob- 
jects with a long-lived disk, this value can be -0.01 times lower 



dJura & Kahandll999 ). It is not clear whether M 1-92 has such 
an extra ordinarily low value or not. However, the calculated disk 
mass is at least unlikely to be significantly higher than the afore- 
mentioned values. With models with small grains (micron size 
or smaller) in the disk, it is hard to explain the estimated CO 
disk mass. Taking into account the possible uncertainty by a fac- 
tor of ~2 in the gas-to-dust mass ratio, we conclude that dust in 
the M 1-92 disk is significantly larger than in the envelope and 
the a max value is expected to be larger than 100.0 /mi. 



3. Discussion 

3.1. Grain size 

For the last 20 years, grain growth in AGB stars and PPNs 
has been studied and several observational evidences has been 
reported. Of these, a particular interest to us lies in a series 
of studies in carbon stars that show silicate dust features such 
as the Red Rectangle, BM Gem, V 778 Cyg, and AC Her. 



In this class of objects, the dual chemistry is explained by 
the presence of a disk, in which oxygen-rich dust is stored 
(IWaters et al.l Il998l: iLlovd Evansl Il990b iBarnbaum etaD I1991L 
and references therein). In fact, some observations have de- 
tected disk-like mo t ions in CO emis sion line observations (e.g. 
iKahane et aDll998b iFong et al.l I2006I) and a s patially resolved 
Keplerian ro tating motion in the Red Rectangle dBujarrabal et al.l 
l2003l |2005|) . If such disks live for a long time, dust parti- 
cles can coagulate by grain-grain collis ion and grow in size 
(lJura et al.ll2000bl:lYamamura et al.ll200Qh . More recently, afore- 
mentioned evidence has been also fo und even in some sin- 
gle chemistry systems. For example, Ide Ruyter et al.l (|2006|) 
found flux excesses in A ^ 2.0 jum in their 51 post- AGB 
stars which are known to have binary companions. These in- 
frared excesses indicate small inner radii of dust structures, 
which are kept closer to the star due to their rotating mo- 
tion of the disk rather than their expanding motion. They con- 
clude that presence of disks are common phenomenon in bi- 
nary post- AGB stars. Evidence for large gr ains has also bee n 
found in some PPN s such as AFGL 268 8 (lJura et al.|[2000ab . 
IRAS 17150-3224 dMeixner et alJ 12002b. IRAS 1 8276-1431 
(ISanchez Contreas et al.l |2007|) IRAS 1947 5+3119 dSahai et al] 
2007|), IRAS 22036+5306 dSahaietal] |2008), and IRAS 



222 72+5435 (lUetaet alJ 1200JJ) including our result of M 1- 
92. iFong et al l (120061) performed a survey observation of gas 
phase kinematics of 38 AGBs to PNs, which include some 
of aforementioned objects. The objects with evidence of grain 
growth, except IRAS 22272+5435 (and IRAS 19475+3119), 
have a disk-like motion or a bipolar appearance in the opti- 
cal and NIR images and their initial masses are expected to 
be higher than those of spherical and elliptic PPNs and PNs 
with e xpanding motions (o utflo w) instead of ro tating motions 
(disk) (ICorradi & Schwarzll 19951: iMeixner et alJl2002h . It is ob- 
vious that grain growth is affected by factors such as disk ge- 
ometries, kinematics of the disk, evolution of the stellar system 
(i.e. duration time of the dust processes), and the stellar masses. 
Previous work has provided qualitative evaluations to determine 
if large grains exist or not and if a disk-like structure exist or 
not to explore the possibility of the disk hypothesis and grain 
growth. Future work should probably be focused on more qual- 
itative, direct discussions on the disk hypothesis. For this pur- 
pose, a more systematic and self consistent analysis are required 
to estimate the physical parameters of the disk and dust. 

With respect to dust, the presence of large dust grains in the 
nebulae is also intriguing. Optical to NIR imaging polarimetry 
offers a great advantage in this issue. In polarimetric images 
of circumstellar dust shells around AGB and post-AGB stars, 
a centro- symmetrically aligned polarization vectors are seen. If 
the dust shell geometry and the chemical composition of the dust 
are known, the degree of polarization would provide a good esti- 
mate of the grain size. Particularly, in an optically thin regime 
where most light from the central star is scattered only once 
by dust in the nebula, the upper limit ca n be provided (e.g. 
iDougados et al.ll990l:IPendleton et aTlll99Qh . It is possible to ob- 
tain a rough estimate even from a single wavelength data, under 
some circumstances, which is an advantage of polarimetric im- 
ages compared to intensity images, visibility data, and SEDs. In 
most cases including our sample, intermediate to high polariza- 
tions (P > 20 %) in the NIR are detected and indicate submicron 
sizes in the shell. The dust in nebulae in large distances from 
the central star is thought to be ejected in the AGB mass-loss 
wind with mainly expanding motions instead of orbiting motions 
around the central stars. Thus, these grains are not expected to 
undergo physical processes such as grain-grain collision and ice 
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formation on the grain core much. Thus, dust in this region is 
likely to preserve the original size distribution, which is deter- 
mined in the dust formation process, as explained by the dust 
formation and stella r wind theories (e.g. iGautschy-Loidl et all 
2004; Hofner 2008). In previous observations, we did not find 
any significant difference in dust sizes or at least in th e polariza- 
tion values in various object classes (e.g. |Jural fl996. UMM07). 
It appears that the initial mass of the stars (or the luminosity) 
and the chemical composition of the stellar system (i.e. oxygen- 
rich or carbon-rich) do not affect the dust particle size in the dust 
formation much. 



3.2. Thermal emission images 

In this section, we present some model images in the MIR and 
FIR. As is expressed in the radiative transfer theory, flux from 
objects in these wavelength regions (thermal emission regime) 
is characterized by a combination of the dust temperature and 
density distributions and the dust absorptive opacity, while op- 
tical and NIR fluxes are governed by scattered light. It is worth 
to take multiple wavelength data into account, which is based 
on different physics, to better constrain the model parameters. 
Unfortunately, observed images of M 1-92 in the MIR and FIR 
are not publicly available at this time. Thus, we use the model 
results to predict what we see in the nebulosity and what param- 
eters can be constrained. 



Figured shows X - 10.0 yum and 100.0 yum images (top and 
bottom panels, respectively) of our selected model results with 
different grain populations of a max = 1.0 yum and 1000.0 /mi in 
the disk (left and right panels, respectively). In the 10.0 yum im- 
ages, the surface brightness of the nebulosity is fainter by a factor 
of <10~ 5 than the peak surface brightness of the central source. 
On the other hand, these are comparable in the 100.0 yum im- 
ages. We find a clear correlation between the model images and 
the dust temperature and density distribution maps (Fig.Q]). The 
appearance of the bright central source in the 10.0 yum images 
is due to the presence of a warm (T^300 K), high-density dust 
region (i.e. inner part of the disk). In fact, such resu lts were also 
reproduced by previously reported modeling (e. g. lLopez et all 
l2000HMeixner et al.ll2002l:IUeta & Meixnerll2003h . As explained 
in detail in these previous papers, in the bipolar case, i.e. opti- 
cally thick case at a given wavelength, the inner part and polar 
region of the dust shell are heated well. Taking into account the 
dust density, MIR photons are mainly emitted from a central, 
compact region, but only a little from an envelope far from the 
central star. We illustrate the effect of dust temperature in the 
100.0 yum images by comparing with the 10.0 yum images. As is 
expected in the radiative transfer theory, FIR photons are emitted 
more efficiently from cold dust (T~30 K) than warm dust. In the 
other words, a longer wavelength image traces generally more 
of the outer region. Therefore, in spite of a high density in the 
disk, the central source is faint in this wavelength. In contrast, 
the superwind shell, where cold dust exists, becomes brighter. In 
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terms of the total flux, i.e. the flux seen in the SED plot, the to- 
tal mass of the geometry components are also important because 
the flux is proportional to the dust mass. The estimated mass of 
the superwind shell is 2.6 Mq, which is much higher than the 
disk mass of 0.15 Mq. The FIR flux is governed much more by 
the flux from the superwind than the disk. The FIR data is use- 
ful for additional constraining of the dust mass of the superwind 
and AGB wind shell; however, it is unrealistic to probe the de- 
tailed structure using the present imagers because of a limitation 
of angular resolutions (a few arcsec). 

The surface brightness of the central source would be use- 
ful to constrain the grain sizes in the disk. The 10.0 jim images 
show a brighter central source in the ci mSiX — 1000.0 yum dust 
model than in the a max = 1.0 /mi dust model. As seen in Fig.Q] 
the 10 jim silicate feature (or opacity) becomes weaker in larger 
grains than in smaller grains. This results in lower extinction in 
larger grains. However, the wavelength dependence of the dust 
opacity becomes similar for large grains. In reality, we do not ex- 
pect sufficient accuracy to determine the grain size with observed 
MIR/FIR images, if the d m2LX value exceeds ~10.0 yum. Another 
possible method, besides the one we proposed in this paper, is 
to use the flux excess in the submillimeter and millimeter wave- 
lengths. In these wavelength ranges, the flux from the superwind 
shell decreases by an order compared to FIR flux, as is expected 
in the SED, but flux from the disk does less if the dust in the disk 
is sufficiently grown in size (the spectral opacity index f3 < 1). In 
further works, one can perform high angular resolution imaging 
in these wavelegnth ranges to better constrain the grain sizes and 
disk masses. 



4. Conclusion 

We performed two-dimensional radiative transfer modeling of 
the dust shells of the bipolar PPN M 1-92. Our modeling ap- 
plied geometries with a disk and bipolar envelope surrounded 
by an AGB wind shell, each of which has different dust char- 
acteristics. The model parameters were constrained by compar- 
ing them with the previously observed SEDs, the intensity and 
polarization images from the HST NICMOS 2 archived data 
(UMM07), and a previous radio observation in the C O emis- 
sion line dBuiarrabaletal.lll998bt lAlcolea et al1l2007l) . With a 
waterlily- shaped hollow envelope, the bright bipolar lobes of M 
1-92 were reproduced. For the dust sizes, we found submicron- 
sized grains (a max = 0.5 fim) in the bipolar lobes. The dust 
size in the disk was constrained with the disk mass, which 
was estimated from the CO emission line data. We obtained 
^max = 1000.0 yum. Although this estimate includes some un- 
certainties in the gas-to-dust mass ratio, the grain size distribu- 
tion, and the geometry of the region considered to be the disk, 
an <2 max ^ 100.0 /mi is hard to explain the disk mass and a sig- 
nificantly large size of > 100.0 jim is expected. We conclude 
that grain growth is likely to occur in the M 1-92 disk. Further 
works including submillimeter and millimeter wavelength imag- 
ing and photometry would provide better interpretations on grain 
growth. The small grains in the bipolar lobes are consistent 
with results of many other AGB and PPNs, suggesting that dust 
formed by AGB mass loss has a similar size, which does not 
depend on the object much. The formation of such submicron- 
sized dust is explained with dust formation and stellar wind theo- 
ries. The presence of large grains in the disk should be explained 
with another mechanism, which is probably the long-lived disk 
hypothesis. The grain growth possibly depends on the disk ge- 
ometry and the stellar temperature. Better understanding will be 



provided if a detailed systematic analysis is made for many other 
samples. 
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